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Abstract 

A  model  to  determine  magnetospheric  electron  fluxes  from  bomb 
neutron  decay  is  presented.  The  source  of  electrons  was  determined 
from  a  neutron  decay  model  dependent  on  a  37  group  neutron  spectrum 
of  a  nuclear  detonation.  Monte-carlo  simulation  of  the  decay  process 
determined  the  electron  energy  and  angular  spectra  as  a  function  of 
neutron  energy.  The  Air  Force  Weapons  Laboratory's  PROMPT  code 
generates  neutron  spectra  from  a  nuclear  detonation  at  an  observer. 
Using  this  neutron  spectrum  and  the  electron  energy  and  angular 
distributions  from  neutron  decay  theory,  an  algorithm  was  developed 
to  calculate  the  electron  flux.  An  example  of  a  1  kt  burst  at  20  km 
shows  that  fluxes  above  normal  background  can  be  observed  as  far 
away  as  400  km.  Higher  bursts  or  bursts  of  larger  magnitude  will 
produce  higher  fluxes  at  this  distance.  It  is  thus  concluded  that 
neutron  decay  alone  can  generate  significant  numbers  of  energetic 
electrons  in  the  magnetosphere.  Their  effect  may  prove  to  be  of 
tactical  significance  in  medium  to  high  altitude  burst  scenarios 
where  fission  debris  is  a  less  significant  source  of  electrons. 
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Chapter  I 
Introduction 

High  altitude  nuclear  bursts  can  inject  large  numbers  of  high 
energy  electrons  and  protons  into  the  magnetosphere.  These  particles 
can  become  trapped  in  the  magnetic  field  and  enhance  the  existing 
radiation  belts.  Satellites  passing  through  these  belts  can  be 
damaged  when  the  radiation  is  severe  enough  to  penetrate  the  shielding 
and  reach  the  large  scale  integrated  circuitry  which  has  a  relatively 
low  radiation  damage  threshold  (Ref  1:1).  Satellites  can  be  damaged 
by  the  lower  energy  electrons  from  surface  charging  in  the  satellite 
dielectrics.  Because  satellites  are  expensive  and  serve  many  impor¬ 
tant  peacetime  and  wartime  missions,  satellite  survivability  and 
vulnerability  concerns  both  private  industry  and  the  government. 

This  study  considers  the  electrons  from  bomb  neutron  decay 
independent  from  all  other  electron  sources.  Generally,  the  most 
significant  source  of  electrons  from  nuclear  bursts  is  the  beta  decay 
of  the  radioactive  fission  debris  nuclei  (Ref  2).  However,  another 
source  of  electrons  is  from  bomb  neutron  decay.  Although  neutron 
decay  creates  significantly  higher  electron  fluxes  than  the  normal 
background,  this  source  is  often  relatively  insignificant  compared 
with  the  electron  fluxes  from  the  fission  debris  (Ref  3:159,  4:4637). 
But  cases  do  exist  where  early  time  electron  fluxes  are  high  and 
cannot  be  explained  by  the  fission  debris  source  (Ref  5:5).  Also 


cases  exist  where  the  atmosphere  severely  attenuates  the  fission 
debris  while  the  neutrons  are  much  less  attenuated.  The  specific 
purpose  of  this  research  was  to  model  neutron  decay  electron  injection 
into  the  terrestrial  magnetosphere  from  nuclear  bursts. 

This  was  done  to  supplement  the  Air  Force  Weapon  Laboratory's 
(AFWL)  capability  to  calculate  electron  environments  in  the 
magnetosphere.  A  prompt  radiation  computer  code  existing  at  AFWL 
calculates  the  neutron,  gamma,  and  x-ray  spectra  from  a  given  nuclear 
detonation  observed  at  an  arbitrary  point  above  the  horizon.  These 
spectra  are  applied  to  satellite  survivability/vulnerability  analysis, 
but  they  lacked  the  electron  spectrum  injected  by  the  decaying 
neutrons.  Such  an  injection  model  was  formulated  and  added  to  AFWL's 
PROMPT  radiation  code. 

The  procedure  to  investigate  this  problem  involves  both  the 
energy  and  the  angular  electron  distributions  from  neutron  decay. 

Using  monte-carlo  techniques  to  model  the  decay  process,  the  electron 
energy  and  angular  distributions  were  determined  as  a  function  of 
neutron  energy.  These  distributions  were  added  as  subroutines  to 
AFWL's  PROMPT  computer  code  to  calculate  the  electron  flux  from 
arbitrary  bomb  neutron  spectra.  To  validate  the  model,  a  test  case 
compared  the  theoretical  results  with  Explorer  IV  satellite  data  for 
the  Teak  event,  data  that  could  not  be  explained  by  fission  debris 
decay.  Other  validation  tests  compared  the  theoretical  results  with 
the  zero  mass  integral  and  spherical  divergence  approximations. 
Finally,  the  dependence  between  mass  integral  and  electron  flux  was 
related  to  the  effect  of  different  altitude  bursts  on  the  electron 
flux  seen  by  a  satellite. 
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Chapter  II 
Background 

23 

Nuclear  bursts  produce  about  10  neutrons  per  kiloton  of  fission 
yield.  These  neutrons  can  interact  with  the  atmosphere,  scattering 
until  they  are  absorbed.  Until  they  are  absorbed,  they  have  a 
finite  probability  of  decaying  into  an  electron,  proton,  and  anti¬ 
neutrino.  Neutrons  have  a  half  life  of  10.6  minutes  (Ref  6),  and 
decay  according  to  equation  1, 

N - ^ -  (1) 

4uR  v  t 
n  n 

where  N  is  the  number  of  neutrons  that  decay  per  unit  volume, 

S  is  the  number  of  neutrons  from  the  source,  R  is  the  distance  from 
the  burst,  v^  is  the  neutron  velocity,  and  is  the  mean  neutron 
lifetime . 

The  atmosphere  readily  absorbs  most  neutrons  and  makes 
neutron  decay  a  minor  source  of  electrons  deep  in  the 
atmosphere.  However,  where  the  atmosphere  is  thin  enough  to  propa¬ 
gate  neutrons,  neutron  decay  can  become  a  significant  contributor  to 
the  electron  environment  (Ref  3:159,  4:4637). 

In  most  electron  trapping  studies,  the  radioactive  fission 


debris  nuclei  are  considered  the  primary  source  of  electrons.  For 
most  applications  this  assumption  is  correct,  because  neutron  decay 
produces  only  one  electron,  while  fission  debris  decay  produces  about 


seven  electrons  (Ref  2).  But  the  actual  relative  significance  will 
be  a  trade-off  between  the  fission  fragment  and  the  neutron  decay 
electron  production. 

In  1962,  Killeen,  Hess,  and  Lingenfelter  estimated  the  electron 

fluxes  caused  by  neutron  decay  (Ref  4).  They  concluded  that  from 

Starfish,  a  1.4  megaton  burst  at  400  km,  the  expected  electron  flux 

7  2 

from  neutron  decay  would  be  less  than  10  electrons/cm  -sec  near 
the  explosion.  This  flux  is  a  small  fraction  of  the  observed  elec¬ 
tron  flux  (Ref  4:3467).  However,  in  reference  to  this  article,  Hess 
later  said  that  this  electron  flux  is  not  negligible,  but  is  con¬ 
cealed  by  the  considerably  larger  flux  of  fission  debris  electrons 
from  Starfish  (Ref  3:159). 

More  recently,  this  author  estimated  the  electron  flux  from 

neutron  decay  as  a  function  of  distance  from  a  burst.  Although  the 

ideal  situation  of  free  space  and  monoenergetic  neutrons  was  examined, 

it  showed  that  under  these  conditions,  neutron  decay  contributes 

significant  electron  fluxes  above  normal  synchronous  background  at 

great  distances  from  a  deep  space  burst  (see  Appendix  A).  The  normal 

6  2 

background  at  synchronous  orbit  is  about  10  electrons/cm  -sec 

(Ref  7),  and  calculations  showed  that  a  one  megaton  burst  in  free 

space  could  generate  this  flux  over  1,000  km  from  the  point  of  burst. 

In  reference  to  these  high  fluxes,  bursts  like  Teak  (75  km)  and 

Orange  (45  km)  should  provide  some  observational  data  since  the  air 

density  is  larger  than  that  for  Starfish  (400  km),  but  low  enough  to 

propagate  neutrons  (Ref  8).  Observation  of  electron  densities  of 
4  2 

4.3  x  10  e“/cm  (Ref  5:10)  were  observed  by  the  Explorer  IV  satellite 
only  18  minutes  after  the  Teak  detonation.  The  spectrum  of  the 


electron  energies  was  not  characteristic  of  fission  debris  (Ref  5:5) , 
and  an  explanation  has  not  been  published  to  this  date.  Some  people 
postulate  that  the  fluxes  could  be  from  protons. 

Explorer  IV  had  3  electron  detectors,  E>580  keV,  E>1  MeV,  and 
E>5  MeV.  These  energy  groups  do  not  provide  an  ideal  spectrum  for 
confirming  or  denying  the  neutron  decay  source,  however,  the 
observed  flux  peak  for  E<1  MeV  is  consistent  with  a  neutron  decay 
electron  spectrum.  AFWL  supported  this  study  to  better  quantify  the 
neutron  decay  source  in  terms  of  magnitude,  energy  spectra,  and 
spatial  distribution  as  a  function  of  burst  yield,  location,  and 
observation  point. 
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Chapter  III 

The  Beta  Spectrum  from  Neutron  Decay 

Electron  flux  in  the  exosphere  is  a  time  dependent  phenomenon 
depending  on  the  trapped  electron  distribution.  The  trapped  elec¬ 
tron  distribution  depends  on  the  electron  motion  in  the  magnetosphere. 
Since  charged  particle  motion  in  a  magnetic  field  depends  on  the 
particle  energy  and  the  angle  between  its  velocity  vector  and  the 
magnetic  field  line,  these  two  parameters  needed  to  be  investigated. 
Since  no  adequate  theoretical  model  to  describe  neutron  beta  decay  was 
available,  data  was  obtained  from  monte-carlo  simulation  of  neutron 
beta  decay.  This  chapter  is  broken  into  three  sections.  The  sec¬ 
tions  present  the  monte-carlo  data  generation,  the  resulting  energy 
distribution,  and  the  resulting  angular  distribution. 

Monte-Carlo  Data  Generation 

No  adequate  theories  exist  to  accurately  describe  the  electron 
energy  and  angular  distributions  from  neutron  decay.  Some  beta 
spectra  of  thermal  neutron  decay  are  available ,  but  theories  do  not 
accurately  predict  the  lower  energy  regions.  Without  an  acceptable 
theory  to  build  a  mathematical  model,  a  monte-carlo  simulation 
of  the  neutron  decay  phenomenon  was  used  to  generate  a  data  base 
for  analysis. 

The  simulation  was  based  on  energy  and  angular  probability 
functions  for  specific  neutron  energies.  A  beta  spectrum  of  thermal 
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neutrons  measured  in  1951  by  Robson  (Ref  9:352)  was  fit  with  a  fourth 
order  polynomial  of  electron  energy  in  keV.  This  curve  was  normal¬ 
ized  and  integrated  so  that  random  numbers  could  be  used  to  specify 
decay  energies.  Electron  energy  groups  of  50  keV  each  from  0  to 
1.15  MeV  were  used  to  accumulate  energy  spectra  histories.  Electron 
spectra  from  37  neutron  energy  groups  were  generated  assuming  that 
the  electrons  were  isotropic  and  that  the  electron  originated  with 
the  same  velocity  as  the  neutron. 

Using  relativistic  velocity  addition  and  randomly  simulated 
isotropic  electron  decay  in  the  center  of  mass  reference  frame,  the 
simulation  calculated  the  electron  spectral  and  angular  data  for 
specified  neutron  energies.  20,000  histories  were  sufficient  to 
calculate  energy  distributions  with  standard  deviations  (/N/N)  on  the 

_3 

order  of  10  for  an  average  group. 

The  angular  distribution  simulation  assumed  isotropic  scatter, 
so  a  two  dimensional  geometry  was  used  to  take  advantage  of  symmetry. 
The  x-axis  was  the  direction  of  the  neutron  velocity  vector  and  the 
y-axis  was  perpendicular  to  the  x-axis.  Realizing  that  the  angular 
distribution  relies  on  the  third  dimension,  the  probability  of  an 
electron  emitting  at  an  angle  6  above  or  below  the  x-axis  is 

P<e>  . 

where  6  is  defined  from  0  to  •&.  This  is  the  distribution  function 
used  to  randomly  generate  6  in  the  center  of  mass  reference  frame. 

Once  theta  was  generated,  the  electron  energy  was  randomly 
generated.  The  combination  of  this  angle  and  energy  was  then  added 


relativist ically  to  the  neutron  velocity  to  determine  the  lab  frame  elec¬ 
tron  energy  and  angle  distributions.  The  data  was  tabulated  in  electron 
energy  groups  of  SO  keV  each  and  in  angular  groups  of  5  degrees  each. 

The  resulting  angular  information  depended  on  neutron  energy  as 
well  as  electron  energy.  Therefore,  for  every  50  keV  electron  energy 
group,  the  angular  information  was  divided  into  38  angular  groups. 

To  account  for  the  increased  number  of  data  groups,  the  iterations 
were  increased  to  500,000  histories  to  get  a  standard  deviation  of  10  . 

The  monte-carlo  simulation  was  validated  by  comparing  two  neu¬ 
tron  energy  extremes  to  the  original  data.  A  zero  energy  neutron  case 
was  tested  and  generated  the  same  distribution  that  was  input.  A 
20  MeV  neutron  case  generated  an  angular  spectrum  skewed  in  the 
forward  direction  and  a  considerably  higher  energy  skewed  electron 
spectrum.  Runs  using  other  sets  of  random  numbers  were  examined  and 
observed  deviations  were  within  expected  random  error. 

Polynomial  curve  fits,  based  on  the  monte-carlo  data,  were  then 
used  to  save  computer  space  and  time  during  parametric  studies. 

Each  set  of  data  for  the  electron  energy  spectrum  was  curve  fit  by 
a  fourth  order  polynomial  equation  in  electron  energy  in  keV.  Each 
set  of  data  for  the  angular  distribution  was  curve  fit  using  a 
fourth  order  polynomial  equation  in  theta  in  radians.  The  following 
sections  describe  each  topic  in  greater  detail. 

Electron  Energy  Distribution  from  Neutron  Decay 

To  determine  the  electron  flux,  the  electron  velocity  must  be 
known.  Electron  velocities  are  calculated  from  their  energies 
relativistically,  so  the  velocity  distribution  is  highly  dependent 
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on  the  energy  distribution.  This  section  presents  the  theory  and 
modelling  of  the  electron  energy  distribution  from  an  arbitrary 


neutron  spectrum. 

The  free  neutron  beta  decay  spectrum  has  been  investigated  for 
many  years,  but  unLil  recently  the  equipment  required  to  observe  the 
electron  spectrum  from  free  neutrons  near  rest  has  not  been  available. 
Golub,  et  al.,  proposed  research  in  this  field  by  using  cryogenic 
bottles  to  contain  free  neutrons  for  observation  (Ref  10:134). 
Unfortunately,  this  data  has  not  been  published  yet,  and  could  not 
be  used  in  this  study.  However,  the  neutron  beta  decay  spectrum 
measured  in  by  Robson  is  reliable  and  is  still  being  used 

today  (Ref  9:3b2,  11:378).  It  was  thus  chosen  for  this  study. 

Figure  1  shows  Robson's  thermal  neutron  decay  beta  spectrum  with  his 
curve  fit  for  the  data.  Since  no  information  was  available  on  the 
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Fig.  1  The  Experimentally  Observed  Neutron 
Beta  Decay  Energy  Spectrum 
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lower  energy  electrons,  this  author  speculated  the  zero  energy 

2 

intercept  and  used  the  .4  moc  (mQ  is  the  electron  rest  mass)  and 

higher  energy  data  points  to  generate  coefficients  for  a  polynomial 

curve  fit.  A  fourth  order  polynomial  fit  was  shown  to  fit  the  data 

with  a  Chi-square  value  of  .2521,  which  is  off  by  less  than  0.5 

percent.  The  unnormalized  polynomial  fit  is  shown  in  equation  2 

2 

with  electron  energy  (in  mQc  )  as  the  independent  variable. 

2  3 

P(E)  *.0268061907  +  5.554304753E  -  9.207097299E  +  5.101880231E 

.9502414486E4  (2) 

Figure  2  shows  the  superposition  of  this  curve  fit  normalized 
with  the  experimental  data.  Robson's  data  shows  the  maximum  electron 
energy  at  844  keV.  Theoretically,  the  maximum  electron  energy 
should  be  785  keV.  The  higher  value  is  within  experimental  error 


in 


0.0  287.5  575.0  862.  S  1150.0 

ENERGY  (KEV) 

Fig.  2  Beta  Energy  Spectrum  from  Neutron  Decay 
Curve  Fit  Superimposed  Over  Data 
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of  the  theoretical  limit,  and,  not  knowing  the  nature  of  the  error, 
the  errors  caused  by  forcing  the  data  to  fit  the  theoretical  limit 
would  be  unpredictable  and  possibly  invalidate  the  entire  spectrum. 
Therefore,  the  higher  energy  electrons  were  included  in  the  model. 

This  beta  spectrum  will  depend  on  the  neutron  energy. 

Assuming  that  the  electron  and  proton,  decay  products,  possess  the 
same  velocity  as  the  neutron  before  it  decays,  the  neutron  velocity 
must  be  added  relativistically  to  the  velocity  each  particle  obtains 
from  the  decay  energy.  Since  we  are  considering  relativistic 
particles,  velocity  addition  dramatically  increases  the  electron 
energy  in  order  to  conserve  momentum.  Hence,  increasing  neutron 
energy  skews  and  squashes  the  electron  energy  spectrum.  Figure  3 
shows  the  difference  between  the  beta  spectra  of  0  and  19.64  MeV 
neutrons.  These  two  curves  represent  the  two  extremes  expected  for 
neutron  energies  from  bomb  neutrons. 

Fourth  order  polynomial  curve  fits  were  calculated  for  each 
neutron  energy  group,  and  then  each  fit  constant  was  found  as  a  function 
of  neutron  energy.  The  general  equation  determining  the  electron  spectra 
from  an  arbitrary  neutron  spectrum  is  shown  in  equation  3.  Appendix  B 
presents  the  polynomial  fit  constants  of  the  beta  spectrum  curve  fits 

N(E  ,E  )  -  bn(E  )  +  b.(E  )E  +  b„(E  )E  *  +  b,(E  )E  3  + 
n  e  On  lne  2ne  3ne 


wv 

4  n  e 


(3) 
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0.0  387.  S  575,0  862.5  1150.0 

ENERGY  (KEV) 

Fig.  3  Normalized  Beta  Energy  Spectra  for 
0  MeV  and  19.64  MeV  Neutrons 

for  various  neutron  energies.  The  neutron  energy  groups  were  chosen 
to  coincide  with  the  important  neutron  energy  groups  for  calculating 
the  neutron  induced  secondary  gammas  since  they  were  already  used  in 
the  PROMPT  code  (Ref  12:7). 

The  curves  fitting  the  fit  constants  for  arbitrary  neutron 
energies  is  also  given  in  Appendix  B.  These  constants  correspond  to 
the  in  equation  3.  Since  the  neutron  energy  groups  are  narrow, 
linear  interpolation  between  the  fit  constants  could  be  used  as  a 
convenient  way  to  calculate  the  coefficients.  Equation  3,  however, 
provides  slightly  more  accurate  results  with  less  computer  space. 

An  example  compares  the  expected  curve  fit  found  with  equation  3  to 
the  linear  interpolation  between  the  12.82  MeV  and  the  13.84  MeV 
neutron  energy  groups.  A  neutron  energy  of  13.00  MeV  is  used  for 
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comparison.  To  determine  the  error,  another  monte-carlo  run  was 
made  for  this  neutron  energy.  The  difference  between  the  equation  3 
result  and  the  monte-carlo  result  is  .0855,  while  the  difference  was 
.12%  for  linear  interpolation. 

Since  the  equations  are  normalized,  the  electron  spectrum  from 
any  type  of  neutron  spectrum  can  now  be  found.  Arbitrary  electron 
energy  groups  can  be  defined  to  acquire  information  on  the  electron 
spectrum  from  bomb  neutron  decay.  These  groups  are  defined  between 
the  energies  and  E^+^.  Integrating  equation  2  determines  the 
expected  number  of  electrons  per  source  neutron  between  the  energies 
and  E^+^.  This  result  can  be  multiplied  by  the  number  of  neutrons 
that  decay  in  that  neutron  energy  group  and  then  summed  over  all 
neutron  groups.  Equation  4  summarizes  this  calculation. 


m 

Z 

j=0 


4 

Z 

£=0 


a‘<Ei)<E«i.f  V 


£+1 


7I+iT 


(4) 


In  equation  4,  Neg.  is  the  number  of  electrons  with  energy  between 
E^  and  E.^.  Nng •  *s  num^er  of  neutrons  with  energy  E^  of  the 
j  th  neutron  energy  group  that  decay,  and  a^  are  the  appropriate 
polynomial  curve  fit  coefficients  found  in  Table  B-I  of  Appendix  B. 


Electron  Angular  Distribution  from  Neutron  Decay 

In  addition  to  the  electron  energy  distribution,  the  electron 
angular  distribution  must  be  known  to  model  the  trapped  electron 
distribution.  This  section  presents  the  theory  and  modelling  of  the 
electron  angular  distribution  from  neutron  decay. 
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Theoretically,  the  electron  energy  limit  imposed  by  the  decay 
energy  is  785  keV.  Any  electron  observed  at  a  higher  energy  must 
originate  from  an  energetic  neutron.  The  fastest  electron  possible 


from  a  decay  must  be  scattered  forward  at  zero  degrees  from  the 
neutron  velocity  vector.  The  slowest  electron  will  be  one  emitted 
opposite  to  the  neutron  velocity  vector.  This  indicates  that  the 
electron  angular  distribution  depends  on  both  the  electron  energy 
and  the  neutron  energy. 

The  final  electron  energy  distribution  from  the  decay  is 
different  for  different  neutron  energies.  This  correlation  was 
derived  in  equations  2  and  3.  Each  spectrum  is  different  because 
the  neutron  adds  velocity  to  the  electron  and  tends  to  skew  the 
electron  angular  distribution  forward,  as  shown  in  figure  4.  How¬ 
ever,  the  electron  angular  distribution  also  depends  on  the  electron 


energy  distribution  because  a  14  MeV  neutron  can  create  a  1  MeV 
electron.  However,  to  conserve  momentum  and  energy,  the  1  MeV 
electron  will  only  be  observed  at  an  angle  very  forward  scattered. 
Modelling  this  information  independent  of  electron  energy  might  be 
useful  to  observe  the  scatter  distribution  of  the  electrons  as 
figure  4  shows.  But,  the  important  information  is  the  relationship 
between  electron  energy  and  its  angular  distribution.  This  relation¬ 
ship  depends  on  both  the  electron  and  the  neutron  velocities. 

Appendix  C  shows  the  angular  distribution  of  0-50  and  500-550  keV 
electrons  for  various  neutron  energies.  It  assumes  isotropic  decay 
in  the  center  of  mass  reference  frame. 

Assuming  that  the  decay  of  neutrons  at  rest  is  isotropic, 
the  same  general  monte-carlo  method  that  generated  the  energy 
spectra  simultaneously  generated  the  angular  information. 

Simultaneous  generation  is  critical  since  the  electron  angular 
distribution  greatly  influences  the  electron  energy  distribution 
for  a  specific  neutron  energy  as  found  in  the  previous  section. 

The  electron  angular  distribution  depends  on  neutron  energy  and 
electron  energy.  However,  figure  4  shows  the  angular  distribution 
for  all  electrons,  independent  of  electron  energy,  per  source 
neutron  of  a  zero  and  a  19.64  MeV  neutron. 

For  electrons  originating  from  0  to  19.64  MeV  neutrons,  the 
electron  spectrum  ranges  from  0  to  1.1  MeV.  This  electron  energy 
range  was  broken  into  50  keV  energy  groups  for  angular  distribution 
modelling.  Fourth  order  curve  fits  were  found  for  each  50  keV  group 
for  12  neutron  energies  to  solve  equation  5,  where  6  is  in  radians. 
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(5) 


P(0)  “  a_(E  )  +  a.(E  )0  +  a  (E  )02  +  a,(E  )03  +  a,(E  )04 
Un  in  i  n  j  n  4  n 


This  data  base  was  very  large,  so  a  7th  order  curve  fit  for  each  a^ 
dependent  on  neutron  energy  was  found.  Equation  6  is  used  to 
calculate  the  a^  using  the  fit  constants  presented  in  subroutine 
ANGLE  of  Appendix  E. 


a. (E  )  =  b_  +  b.E  +  b.E  2  +  b,E  3  +  b.E  4  +  b.E  5  +  b,E  £ 
i  n  0  In  2n  3n  4n  5n  6n 


b.E  ' 
7  n 


(6) 


where  the  E  are  in  keV. 
n 

From  this  data  base,  the  distribution  of  electron  angles  for 
electrons  in  specific  50  keV  groups  can  be  found.  The  a^  fit 
coefficients  for  a  specific  50  keV  energy  group  can  be  found  for 
each  neutron  energy  group.  These  coefficients  are  used  to  calculate 
the  electron  angular  distribution  by  integrating  equation  5  and 
defining  angular  bins  to  any  specified  detail. 

Since  the  electron's  angle  with  the  magnetic  field  line  is 
important  to  calculate  the  particle  motion,  a  correlation  between 
the  decay  angle  and  the  magnetic  field  line  is  required.  Appendix  D 
derives  the  relation  between  a  neutron  incident  on  a  magnetic  field 
line  and  the  electron  density  as  a  function  of  an  angle  a  from  the 
magnetic  field  line.  The  procedure  outlined  in  Appendix  D  can  be 
used  to  numerically  evaluate  the  number  of  electrons  between  a  and 
a  +  da  from  a  field  line.  This  technique  requires  transforming 
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the  coordinates  of  equation  5  into  coordinates  defined  around  the 
magnetic  field  line. 
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Chapter  IV 

Distribution  Modelling 


The  Air  Force  Weapons  Laboratory  computer  code  PROMPT  calculates 
neutron,  gamma,  and  x-ray  spectra  from  nuclear  detonations  for  points 
in  space  from  the  prompt  radiation.  Since  electrons  create  hazards 
in  space,  the  capability  to  calculate  the  electron  spectrum  and  flux 
at  a  point  in  space  from  neutron  decay  needed  to  be  added.  This 
chapter  is  devoted  to  how  the  theory  developed  in  the  previous 
chapters  was  added  to  PROMPT.  The  subroutines  are  provided  in 
Appendix  E  with  a  user's  manual. 

Electron  Energy  Spectrum 

The  electron  energy  spectrum  developed  early  in  Chapter  III 
depends  on  the  bomb  neutron  energy  spectrum.  Since  the  PROMPT  code 
already  calculates  the  neutron  spectrum  at  a  point,  this  information 
provides  the  electron  density  and  the  electron  energy  distribution 
from  each  neutron  energy  group.  Taking  the  electron  spectrum  curve 
fit  constants  as  determined  by  the  neutron  energy  applied  to  the 
constants  of  Table  B-II  in  Appendix  B,  the  electron  energy  spectrum 
for  each  of  the  38  neutron  energy  groups  is  found.  Using  50  keV 
electron  energy  groups  for  spectral  shaping,  23  electron  groups 
between  0  and  1.15  MeV  contain  the  electron  density  distribution 
from  neutron  decay. 

■« 

•r* 
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Applying  equation  1  to  the  neutron  density,  and  taking  into 
account  the  time  dilation  of  the  half  life  of  the  higher  energy 
neutrons,  the  electron  density  from  neutron  decay  is  found.  This 
calculation  uses  the  lower  neutron  energy  from  each  neutron  energy 
group.  Taking  each  neutron  energy  group  separately,  the  electron 
energy  spectrum  is  integrated  and  the  number  of  electrons  in  each 
50  keV  electron  energy  group  is  found  by  multiplying  the  number  of 
decayed  neutrons  by  the  expected  number  of  electrons  in  each  electron 
energy  group  per  source  neutron.  The  total  number  of  electrons  in 
each  of  these  energy  groups  is  found  by  summing  the  contributions 
from  all  neutron  energy  groups.  The  resulting  distribution  is  the 
desired  total  electron  energy  spectrum. 


Electron  Angular  Distribution 

The  electron  angular  distribution  developed  in  the  later  part 
of  Chapter  III  shows  that  the  angular  distribution  depends  on  both 
the  neutron  and  the  electron  energies.  For  a  point  in  space,  this 
subroutine  was  developed  to  calculate  the  angular  distribution  with 
respect  to  the  forward  scattered  angle  of  0  radians.  The  method  to 
calculate  this  distribution  begins  by  calculating  the  curve  fit 
coefficients,  the  a^  in  equation  7,  of  the  electron  angular  distribu¬ 
tion  for  each  neutron  energy  group.  Equation  7  requires  the  neutron 
group  energy  E^  and  the  polynomial  fit  coefficients  b^  of  the  i1-*1 
electron  energy  group.  The  resulting  angular  distribution  of  elec¬ 
trons  decaying  from  a  specific  energy  neutron  are  modelled  with 
fourth  order  curve  fits.  For  a  specific  electron  energy,  this 

distribution  determines  the  N_  (0)  in  equation  8.  This  distribution 

e 
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is  applied  to  each  of  the  electron  energy  groups  previously  mentioned, 
except  that  this  distribution  depends  on  the  0-50,  50-100,.., 

1100-1150  keV  electron  energy  groups  specifically  because  the  curve 
fit  constants  were  derived  dependent  on  these  groups. 


a. (E  )  =  I  b..E  J 
1  n  j=0  J1  n 


(7) 


Ne  (6) 

e,n 


4 

I  a.  e1 
_ a  in 


(8) 


Each  of  these  fourth  order  equations  are  integrated  to  determine 
the  distribution  of  electrons  between  specific  angular  groups.  The 
total  number  of  electrons  are  then  multiplied  by  the  angular  distri¬ 
bution  function  for  electrons  just  determined.  This  method  is 
repeated  for  each  neutron  energy  group  until  all  the  electrons  from 
each  neutron  group  are  distributed  into  the  angular  bins.  Equation  9 
presents  the  calculation  of  the  total  angular  distribution  of  elec¬ 
trons  with  respect  to  the  angle  from  the  forward  neutron  velocity. 
N(0)  is  the  distribution  of  electrons  at  an  angle  0,  Ng  is  the  total 
number  of  electrons,  En^  is  the  neutron  energy  of  the  i^h  group,  and 
the  other  parameters  are  consistent  with  equations  7  and  8. 


N(0)  -  N 

e 


37 

I 

i=l 


(9) 


The  angular  information  is  important  to  determine  the  electron 
life  span  in  the  magnetosphere.  Since  the  flux  depends  on  electron 
energy,  by  virtue  of  the  velocity,  the  flux  will  vary  with  time  as 


the  electrons  leave  the  trapped  region.  Therefore,  the  angular 
distribution  of  electrons  with  specific  energies  relative  to  the 
magnetic  field  must  be  known.  Because  of  this  problem,  as  well  as 
the  problem  that  the  angular  distribution  depends  on  both  the  elec¬ 
tron  and  neutron  energy,  the  total  angular  distribution  is  of  little 
applicable  use,  but  the  subset  of  distributions  are  critical  for 
flux,  time  dependent  models.  Appendix  E  addresses  the  numerical 
solution  to  this  part  of  the  problem. 

One  should  note  that  PROMPT  provides  the  neutron  spectra  at  a 
point  in  space,  but  it  does  not  provide  any  information  concerning 
the  angular  nature  of  the  neutrons.  Due  to  scatter  in  the  atmosphere, 
some  neutrons  may  not  arrive  radially  from  the  burst.  In  a  homoge¬ 
neous  atmosphere,  however,  one  would  expect  the  average  neutron  to 
arrive  radially  from  the  burst.  This  means  that  the  results  of  this 
study  would  hold  true  in  an  average  sense.  An  inhomogeneous  atmos¬ 
phere  would  be  impossible  to  handle  with  PROMPT  because  the  neutron 
angular  distributions  require  significantly  more  sophisticated 
techniques.  Assuming  radial  neutrons  serves  as  a  source  of  error  in 
the  electron  angular  distributions,  but,  given  the  correct  neutron 
angular  information,  the  theory  and  method  could  be  simply  extended. 
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Chapter  V 

Flux  Calculations  in  the  Magnetosphere 

Foundations 

Since  the  discovery  of  trapped  radiation  in  the  magnetosphere, 
many  people  have  wondered  about  the  underlying  physical  properties 
involved  in  this  phenomenon.  Experiments  designed  to  study  particle 
motions  in  the  magnetosphere  used  charged  particle  accelerators  as 
well  as  nuclear  weapons  to  inject  large  numbers  of  charged  particles 
into  the  magnetosphere.  Probably  the  world's  most  noted  experiment 
of  this  sort,  Starfish,  injected  numerous  charged  particles  that 
became  trapped  in  the  magnetosphere  for  years. 

Modelling  the  injection  process  for  bursts  like  Starfish  begins 
with  magnetohydrodynamic  calculations  of  the  plasma  interaction  with 
the  magnetic  field.  From  this  interaction,  the  fission  debris 
distribution  is  found,  and  electron  injection  from  this  debris  is 
calculated.  Electrons  move  between  mirror  points  and  drift  around 
the  earth,  forming  electron  belts.  Some  belts  are  more  stable  than 
others  because  decay  processes  like  atmospheric  absorption,  colli¬ 
sions,  plasma  instabilities,  and  many  other  complicated  mechanisms, 
decrease  the  electron  density  in  these  belts  with  time  (Ref  2). 

The  purpose  here,  is  to  calculate  prompt  electron  fluxes  from 
neutron  decay  electrons  injected  into  the  atmosphere  or  magnetosphere. 
The  source  of  decaying  neutrons  is  prompt  neutrons  from  the  fission 
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and  fusion  processes.  Instead  of  plasma  interactions,  particle 
transport  methods  propagate  the  neutrons  through  real  atmosphere  to 
find  the  neutron  spectrum  at  a  point  in  space. 

The  neutron  transport  method  for  this  investigation  uses  mass 
integral  scaling  fits  of  two  dimensional  anisotropic  neutron  trans¬ 
port  solutions.  These  equations  are  part  of  the  PROMPT  computer 
code  owned  and  maintained  by  the  Air  Force  Weapons  Laboratory  at 
Kirtland  Air  Force  Base  in  New  Mexico.  For  medium  to  high  altitude 
bursts,  the  neutron  transport  calculations  are  in  error  by  less  than 
a  factor  of  two  (Ref  13).  PROMPT  calculates  neutron,  gamma,  and 
x-ray  spectra  at  a  point  in  space.  The  additional  capability  to 
calculate  electron  spectra  from  the  decaying  neutrons  was  added 
during  this  study  to  round  out  the  code  capabilities. 

Electron  Flux  Parameters 

Electron  flux  depends  on  the  number  of  electrons  passing  through 
a  region  of  space  per  unit  time.  The  essential  information  for  flux 
calculations  can  be  found  in  the  electron  spectrum  and  the  electron 
density.  However,  in  the  magnetosphere,  the  electron  interaction 
with  the  magnetic  field  tends  to  reduce  the  electron  density  with 
time  and  to  change  the  spectrum.  For  the  purpose  of  this  paper,  the 
object  is  to  calculate  the  prompt  electron  flux.  Time  dependencies 
were  not  calculated  since  only  the  source  was  being  investigated. 

Techniques  have  been  published  to  calculate  fluxes  along  field 
lines,  (Ref  14)  but  the  electron  spectrum  and  density  must  be  known 
for  each  point  along  the  field  line.  Since  the  foundation  of  the 
prompt  calculation  is  the  mass  integral,  a  relation  between  the  mass 
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integral  and  the  electron  flux  was  found.  The  flux  is  directly 
proportional  to  the  number  of  neutrons,  and  if  one  assumes  the 
spectrum  does  not  change  with  yield,  the  relation  between  flux  and 
yield  is  linear. 

Using  the  PROMPT  code  with  an  unclassified  thermonuclear  spec¬ 
trum,  the  electron  fluxes  as  a  function  of  mass  integral  for  1 
kiloton  bursts  were  calculated  and  are  presented  in  Table  I.  The 
method  used  to  calculate  the  values  in  Table  I  first  chose  a  1  kt 
burst  at  the  altitude  in  the  left  column.  The  electron  flux  was  then 
calculated  at  a  point  400  km  above  the  burst  point.  PROMPT  calcu¬ 
lated  the  mass  integral  in  column  two  and  the  DECAY  subroutine, 
developed  in  this  study,  calculated  the  fluxes  in  column  three. 

Figure  5  plots  the  electron  flux  dependence  on  mass  integral  from  Table  I. 


Table  I 

Electron  Flux  as  a  Function  of  Mass  Integral 


Burst  Altitude 

Mass  Integral 

Electron  Fli 

km 

gm/ cm^ 

e~/cm^-sec 

2 

8.201E+02 

3.935E-09 

5 

5.530E+02 

1.366E-04 

10 

2.686E+02 

1.025E+01 

12 

1.978E+02 

1.957E+02 

15 

1.237E+02 

4.392E+03 

20 

5.590E+01 

5.502E+04 

25 

2.556E+01 

1.117E+05 

27 

1.885E+01 

1.171E+05 

30 

1.208E+01 

1.124E+05 

40 

3.049E+00 

7.989E+04 

50 

8.584E-01 

6.721E+04 

60 

2.380E-01 

6.536E+04 

75 

2.652E-02 

6.703E+04 

100 

3.598E-04 

6.703E+04 

— 

0.0 

5.823E+04 

I 


Table  I  shows  Chat  the  flux  does  not  necessarily  monotonically 

increase  as  the  mass  integral  decreases.  As  the  burst  altitude 

increases,  the  mass  integral  decreases,  but  the  neutrons  have  a 

greater  tendency  to  reflect  off  the  denser  lower  atmosphere.  This 

phenomenon  is  a  combination  of  albedo  and  build  up  and  is  caused  by 

neutron  scatter  in  the  atmosphere.  The  scattered  neutrons  increase 

the  neutron,  and,  thus,  the  electron  density  or  flux.  This  albedo 

effect  has  its  highest  influence  below  40  km.  Outside  this  range, 

the  albedo  neutrons  are  either  not  produced  or  are  negligible  because 

2 

of  spherical  divergence.  With  no  atmosphere,  0  gm/cm  ,  no  build  up 
or  albedo  is  observed.  Hence,  the  last  flux  value  in  Table  I  is  less 
than  the  other  values. 

Comparing  the  fluxes  of  Table  1  to  the  expected  fluxes  calculated 

for  a  vacuum  in  Appendix  A,  Table  I  fluxes  are  higher  than  Appendix  A 

2 

fluxes  until  the  mass  integral  goes  over  100  gm/cm  .  The  discrepancy 

here  lies  in  the  albedo  neutrons,  which  increase  the  electron  density, 

and  a  higher  average  electron  energy  than  that  used  in  Appendix  A. 

These  two  factors  combined  are  enough  to  increase  the  observed  flux  by 

the  appropriate  amount  because  albedo  alone  can  increase  the  neutron 

density  by  a  factor  of  2  or  more  depending  on  the  mass  integral. 

Another  way  to  examine  the  mass  integral  effect  on  electron  flux 

is  to  take  an  observer  at  a  constant  altitude  and  change  the  height 

of  the  burst  directly  below  it.  For  this  example  a  500  km  observer 

was  chosen.  The  burst  was  1  kt  directly  below  the  target  at  various 

2 

altitudes.  Table  II  presents  the  mass  integral  in  gm/cm  and  the 

2 

electron  flux  in  electrons/cm  -sec  for  each  burst  altitude. 

The  information  presented  in  Table  II  is  useful  to  show  the 
general  tendency  that  very  low  altitude  bursts  make  little 
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Electron  Flux  Dependence  on  Mass  Integral 
for  a  Constant  400  km  Distance 


Table  II 

Electron  Flux  as  a  Function  of  Mass  Integral 
for  a  500  km  Observer 


Burst  Altitude 
km 

2 

5 

10 

12 

15 

20 

25 

27 

30 

40 

50 

60 

75 

100 


Mass  Integral 
gm/cm^ 

8.201E+02 
5.530E+02 
2.686E+02 
1.978E+02 
1.237E+02 
5.590E+01 
2.556E+01 
1.885E+01 
1 .208E+01 
3.049E+00 
8.584E-01 
2.380E-01 
2.652E-02 
3.598E-04 


Electron  Flux 
e”/cm^-sec 

2.539E-09 

8.920E-05 

6.832E+00 

1.315E+02 

2.987E+03 

3.821E+04 

7.919E+04 

8.375E+04 

8.140E+04 

6.041E+04 

5.311E+04 

5.402E+04 

5.938E+04 

6.703E+04 
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contribution  to  the  electron  environment  from  neutron  decay.  How¬ 
ever,  it  also  shows  that  bursts  have  increasing  effect  on  the  elec¬ 
tron  environment  up  to  about  30  km  due  to  albedo  and  build  up. 

Bursts  above  30  km  have  decreasing  effect  as  the  albedo  and  build  up 
factors  have  smaller  influence  since  the  mass  integral  decreases. 
When  the  mass  integral  becomes  small  enough  that  neutron  attenuation 
is  negligible,  between  50-60  km,  the  electron  flux  increases  since 


One  other  parameter  essential  for  magnetospheric  flux  calcula¬ 
tions  is  the  electron  angular  distribution  with  respect  to  the 
magnetic  field  line.  Electron  density  in  the  magnetosphere  decays 
with  time  as  the  electrons  trapped  in  the  magnetic  field  interact 
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with  other  particles  and  fields.  The  loss  cone  is  one  of  the  more 
important  loss  mechanisms  for  the  newly  formed  electron  belts.  The 
loss  cone  can  be  defined  as  the  minimum  angle  at  which  electrons  can 
intersect  the  magnetic  field  line  and  mirror.  Electrons  with  angles 
less  than  the  loss  cone  angle  will  not  be  trapped.  Instead,  they  are 
absorbed  by  the  atmosphere.  This  loss  mechanism  rapidly  decreases 
the  electron  density  along  a  field  line  and  decreases  the  observed 
electron  flux.  Although  calculating  the  actual  flux  after  the  elec¬ 
trons  have  dispersed  along  the  field  line  is  beyond  the  scope  of  this 
thesis,  the  angular  distribution  of  the  electrons  with  respect  to  the 
field  line  is  calculated  and  the  underlying  method  is  presented  in 
Appendix  E. 

Each  of  these  parameters  show  why  the  observed  flux  should  be 
lower  than  the  calculated  flux  soon  after  the  prompt  radiation  is 
created.  The  electron  density  is  greatly  dependent  on  slant  range 
due  to  divergence  and  neutron  attenuation.  Other  places  along  a 
magnetic  field  line  could  be  considerably  closer  to  the  burst  than 
where  the  satellite  must  pass  the  field  line.  These  considerations 
show  that  the  prompt  interaction  with  the  entire  field  line  must  be 
calculated  to  determine  the  true  electron  density  in  the  region  that 
the  satellite  passes.  However,  when  the  satellite  observes  the 
prompt  radiation,  the  electron  flux  from  the  prompt  neutrons  should 
be  observed  as  the  theory  predicts. 

Validation  of  the  Theory 

Because  data  on  prompt  radiation  observations  are  unavailable, 
checking  theory  with  observation  is  not  possible.  However,  early 


time  information  was  available  from  the  Teak  event  when  Explorer  IV 

4  2 

reported  fluxes  of  10  (e"/cm  -sec)  for  electron  energies  greater 
than  580  keV  18  minutes  after  the  burst. 

These  fluxes  were  reported  when  the  satellite  was  over  Singapore. 
Unfortunately,  the  position  of  the  satellite  is  not  well-known. 

Some  references  indicate  that  Explorer  IV  was  between  1400  and  1600 
kilometers  in  altitude  at  99°E.  The  Teak  event  occurred  at  17°N 
191°E  and  an  altitude  of  80  km.  The  satellite  could  not  observe  any 
prompt  radiation  because  it  was  below  the  horizon.  The  magnetic 
field  lines  near  the  satellite  never  came  above  the  horizon  either 
(Ref  15).  Thus,  the  electron  fluxes  observed  by  Explorer  IV  must 
remain  an  open  question. 

No  other  early  time  information  was  available  for  any  other 
burst,  so  observational  support  of  the  results  of  this  study  does 
not  exist.  Unless  more  accurate  information  on  the  position  of 

Explorer  IV,  18  minutes  after  Teak,  becomes  available  no  verification 

t 

4  2 

can  be  made.  The  fluxes  of  3  x  10  e"/cm  -sec  could  be  observed  a 
4 

maximum  of  4  x  10  km  from  the  burst. 
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Chapter  VI 

Conclusions  and  Recommendations 


The  purpose  of  this  research  was  to  model  the  neutron  decay 
influence  on  the  electron  space  environment.  During  this  investiga¬ 
tion,  the  neutron  decay  phenomenon  was  broken  into  the  resulting 
electron  energy  and  angular  distribution  as  a  function  of  neutron 
energy.  Assuming  that  the  decay  is  isotropic  in  the  center  of  mass 
frame  and  that  the  neutron  velocity  added  to  the  product  velocities 
from  the  decay  process,  it  was  found  that  the  neutron  energy  greatly 
influenced  both  the  electron  energy  and  the  angular  distributions. 

By  theory,  the  electron  energy  limit  is  only  785  keV,  but,  by  adding 
the  neutron  velocity  relativistically ,  electron  energies  above  1.1  MeV 
could  be  observed.  These  higher  energy  electrons  were  forward 
scattered  as  were  the  rest  of  the  electrons  from  the  higher  energy 
neutrons.  The  existence  of  the  higher  energy  electrons  added 
greatly  to  the  expected  fluxes.  The  actual  damage  these  electrons 
cause  depends  on  the  materials  and  the  shielding  they  interact  with, 
however,  fluxes  greater  than  what  one  would  expect  from  spherical 
divergence  in  free  space  were  calculated  because  of  neutron  scatter 
off  the  earth's  denser  atmosphere.  The  resulting  fluxes  were  higher 
than  normal  background  radiation  in  the  electron  belts  and,  when 
fission  debris  is  not  a  major  source  of  electrons,  neutron  decay 
can  serve  as  a  significant  source  by  itself. 
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By  .determining  that  the  flux  observed  400  km  from  a  1  kt  burst 
5  2 

yields  fluxes  near  10  electrons/cm  -sec,  and  that  the  flux  is 
directly  proportional  to  yield,  any  burst  above  20  km  becomes  a 
significant  source  of  electrons  hundreds  of  kilometers  from  the 
burst  point.  Even  for  a  constant  altitude  observer  at  500  km,  vary¬ 
ing  the  burst  altitude  showed  a  maximum  electron  flux  near  30  km. 

The  albedo  and  build  up  factors  cause  the  different  fluxes  due  to 
the  amount  of  neutrons  they  scatter  towards  the  observer.  Since  the 
electron  flux  is  directly  proportional  to  the  number  of  neutrons,  for 
a  detonation  with  a  similar  neutron  spectrum,  the  electron  flux 
increases  linearly  with  yield.  Cases  were  looked  at  to  explain  the 
Explorer  IV  data,  but  uncertainty  in  the  position  of  the  satellite 
precluded  any  sound  verification  of  the  theory.  No  other  observa¬ 
tional  evidence  was  available  to  compare  with  results  of  this  study. 

The  most  probable  regions  of  errors  would  be  in  the  transport 
and  the  decay  spectra.  Due  to  the  significant  influence  of  the 
neutron  spectrum  on  the  resulting  electron  flux,  the  largest  source 
of  error  would  be  the  neutron  transport  calculation.  The  error  in 
the  neutron  spectrum  for  bursts  above  20  km  is  less  than  a  factor  of 
two,  while  lower  altitude  bursts  can  generate  errors  up  to  one  order 
of  magnitude.  The  major  assumption  underlying  the  theory  of  neutron 
decay  was  the  isotropic  decay  scheme.  If  the  neutron  does  not  decay 
isotropically,  the  electron  energy  and  angular  spectra  would  both  be 
modelled  incorrectly.  The  error  caused  by  this  assumption  is 
probably  negligible. 
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This  thesis  used  purely  analytical  methods  to  model  the  electron 
spectra  as  functions  of  neutron  energy.  The  probability  distributions 
were  complicated  functions,  yet  the  energy  distribution  was  observed 
to  be  similar  to  a  Planckian  distribution  that  flattened  and  skewed 
with  increasing  neutron  energy.  Using  the  neutron  energy  as  a 
temperature  may  simplify  the  modelling  process  and  improve  its 
accuracy.  Presently,  the  only  data  available  on  the  electron  spec¬ 
trum  comes  from  a  1951  experiment.  Theofies  do  not  exist  to  explain 
the  low  energy  electron  spectrum  and  could  be  investigated  using  a 
Planckian  distribution  scheme. 

Until  more  information  is  available  from  satellite  data, 
including  satellite  positions,  and  from  theoretical  development  of 
the  distribution  functions,  no  validations  can  be  made  to  improve 
this  procedure.  Better  approximations  of  the  electron  fluxes  from 
neutron  decay  can  be  made  by  improving  the  transport  calculations, 
the  spectral  information,  or  even  the  curve  fit  solutions,  but  the 
expected  error  is  less  than  a  factor  of  4.  The  necessary  informa¬ 
tion  to  provide  an  accurate  validation  of  this  theory  is  accurate 
early  time  satellite  data  of  a  high  altitude  nuclear  detonation. 

The  general  conclusions  of  this  study  are: 

1.  the  electron  energy  distribution  in  the  lab  frame  is 
greatly  affected  by  the  neutron  energy. 

2.  the  electron  angular  distribution  in  the  lab  frame 
depends  on  the  electron  energy  as  well  as  the  original 
neutron  energy. 

3.  the  electron  energies  from  neutron  decay  are  generally 
in  the  region  influencing  satellite  charging  phenomena. 

4.  the  electron  flux  from  this  source  can  be  significantly 
higher  than  background. 
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Uich  more  research,  the  complete  impact  of  this  source  on  satellite 
survivability/vulnerability  can  be  found.  The  areas  that  need  to  be 
investigated  further  are: 


1.  accurate  satellite  data  including  satellite  position 
of  early  time  electron  flux  observations  of  any  high 
altitude  test, 

2.  investigate  using  the  Planckian  distribution  to  model 
the  electron  energy  spectrum  from  neutron  decay, 

3.  investigate  the  impact  of  this  source  over  an  entire 
magnetic  field  line,  and 

4.  incorporate  the  effect  of  an  incident  neutron  angular 
distribution  on  the  magnetic  field  line  flux. 

Once  these  areas  have  been  studied,  a  complete  analysis  of  the 

problem  should  be  realized. 
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Appendix  A 

Electron  Flux  from  Decaying  Neutrons 

This  appendix  addresses  an  ideal  calculation  of  the 
number  of  neutrons  decaying  at  a  point  in  space  and  the 
electron  flux  at  that  point.  Assuming  no  attenuation  of  the  neutrons 
from  the  burst  and  no  build  up  factor  or  albedo  factor,  the  neutron 
density  that  decays  satisfies  the  equation, 

N - 1 -  (A-l) 

4itR  v  t 
n 

where  N  is  the  number  of  neutrons  per  unit  volume  a  distance  R  from 

a  burst  that  decay,  S  is  the  initial  number  of  source  neutrons,  vr 

is  the  neutron  velocity,  and  t  is  the  neutron  mean  lifetime. 

Additionally,  all  neutrons  are  assumed  to  have  the  same  energy. 

.  23 

For  a  1  kiloton  bomb,  approximately  10  neutrons  are  released 

with  a  mean  energy  of  1.2  MeV.  This  energy  corresponds  to  an  average 

9 

neutron  velocity  of  1.51  x  10  cm/sec.  Since  the  neutron  half  life 
is  10.6  minutes  (Ref  6),  T  is  919  seconds. 

To  determine  the  electron  flux,  equation  A-l  must  be  multiplied 
by  the  electron  velocity.  Using  an  average  electron  energy  of 
600  keV,  the  average  electron  velocity  is  2.66  x  10*^  cm/sec. 

Table  A-l  tabulates  the  expected  electron  flux  as  a  function  of 
distance  from  a  1  kt  burst.  The  fluxes  calculated  for  Table  A-I 
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Table  A-I 


Neutron  Density  and  Electron  Flux 
from  a  1  kt  burst 


R(km) 

N  (e“/cm^) 

2 

Flux  (e"/cm  -sec) 

.1 

5.718E+01 

1.523E+12 

1. 

5.718E-01 

1.523E+10 

10. 

5.718E-03 

1.523E+08 

100. 

5.718E-05 

1.523E+06 

1000. 

5.718E-07 

1.523E+04 

10000. 

5.718E-09 

1.523E+02 

are  highly  spectral  dependent.  Therefore,  the  accuracy  is  limited 
because  of  the  error  from  assuming  all  neutrons  have  the  same  energy. 
Since  the  mean  electron  energy  from  neutron  decay  is  a  little  below 
600  keV,  using  600  keV  as  an  average  energy  keeps  the  error  within 
an  order  of  magnitude.  Therefore,  these  estimates  can  be  used  to 
approximate  the  general  flux  one  might  expect  from  this  situation. 


Appendix  B 

Decay  Electron  Energy  Spectra  Curve  Fitting 

This  section  presents  fit  coefficients  for  the  electron  energy 
spectrum.  Least  squares  polynomial  curve  fits  with  equal  weighting 
of  all  data  points  established  the  equations  discussed  in  this 
section.  Table  B-I  lists  the  fourth  order  fit  coefficients  to  the 
equation: 

P<Ee)  -  a0  ♦  atEe  .  a2Ee2  .  a^3  ♦  a^4  (B-l) 

where  E  is  the  electron  energy  in  keV.  P(E  )  is  the  probability  of 
6  6 

observing  that  energy  electron  and  the  a^  are  the  fit  coefficients 
dependent  on  neutron  energy.  Table  B-II  presents  the  second  order 
fit  coefficients  to  the  a^  in  the  above  equation.  These  coefficients, 
labeled  b^,  satisfy  the  equation: 

a.(E  )  =  bA  +  b,  E  +  b.  E  2  (B-2) 

i  n  0.  1.  n  2.  n 

ii  l 

where  E^  is  the  neutron  energy  in  keV  between  0  and  1.964E+04  keV. 

Special  attention  must  be  made  when  using  these  fits.  They  are 
only  good  for  interpolating  between  0  and  19.64  MeV.  Care  must  be 
taken  to  ensure  that  the  energy  bounds  in  both  equations  are  not 


exceeded. 


E  (MeV) 
n 

1.964E+04 

1.691E+04 

1.492E+04 

1.419E+04 

1.384E+04 

1.282E+04 

1.221E+04 

1.105E+04 

1.000E+04 

9.048E+03 

8.187E+03 

7.408E+03 

6.376E+03 

4.966E+03 

4 . 724E+03 

4.066E+03 

3.012E+03 

2.385E+03 

2.307E+03 

1.827E+03 

1 . 108E+03 

5.502E+02 

]  .576E+02 

1.111E+02 

5.248E+02 

2.479E+02 

2.188E+01 

1.033E+01 

3.355E+00 

1.234E+00 

5.830E-01 

1.013E-01 

2.902E-02 

1.068E-02 

3.059E-03 

1.125E-03 

4.140E-04 

1.000E-08 


Appendix  B 
Table  B-I 


)lynomial  Fit  Constants  of  the  Electron  Spectrum 
for  Different  Neutron  Energies 


a 


4 


7 . 188E-04 
7.092E-04 
7.028E-04 
6.995E-04 
6.995E-04 
6.925E-04 
6.876E-04 
6.842E-04 
6.811E-04 
6.747E-04 
6.725E-04 
6.715E-04 
6.550E-04 
6.630E-04 
6.603E-04 
6.532E-04 
6.493E-04 
6.424E-04 
6.426E-04 
6.428E-04 
6.314E-04 
6.239E-04 
6.137E-04 
6.096E-04 
6.083E-04 
6.077E-04 
6.079E-04 
6.073E-04 
6.055E-04 
6.051E-04 
6.052E-04 
6.051E-04 
6.048E-04 
6 .048E-04 
•050E-04 
.051E-04 
.051E-04 
.050E-04 


1.525E-05 

1.540E-05 

1.562E-05 

1.571E-05 

1.573E-05 

1.589E-05 

1.600E-05 

1.612E-05 

1.624E-05 

1.639E-05 

1.649E-05 

1.655E-05 

1.668E-05 

1.683E-05 

1.688E-05 

1.704E-05 

1.715E-05 

1.730E-05 

1.729E-05 

1.732E-05 

1.753E-05 

1.769E-05 

1.788E-05 

1.796E-05 

1.799E-05 

1.802E-05 

1.802E-05 

1.802E-05 

1.803E-05 

1.807E-05 

1.806E-05 

1.807E-05 

1.807E-05 

1.807E-05 

1 .807E-05 

1.807E-05 

1.807E-05 

1.807E-05 


-5.294E-08 

-5.342E-08 

-5.420E-08 

-5.450E-08 

-5.458E-08 

-5.508E-08 

-5.546E-08 

-5.589E-08 

-5.634E-08 

-5.683E-08 

-5.719E-08 

-5.742E-08 

-5.779E-08 

-5.835E-08 

-5.853E-08 

-5.911E-08 

-5.945E-08 

-5.997E-08 

-5.994E-08 

-6.003E-08 

-6.072E-08 

-6.131E-08 

-6.199E-08 

-6.228E-08 

-6.246E-08 

-6.258E-08 

-6.257E-08 

-6.257E-08 

-6.268E-08 

-6.275E-08 

-6.273E-08 

-6.274E-08 

-6.276E-08 

-6.276E-08 

-6.276E-08 

-6.276E-08 

-6.275E-08 

-6.276E-08 


5.779E-11 

5.822E-11 

5.908E-11 

5.941E-11 

5.949E-11 

6.002E-11 

6.047E-11 

6.093E-11 

6.147E-11 

6.200E-11 

6.241E-11 

6.264E-11 

6.299E-11 

6.364E-11 

6.384E-11 

6.453E-11 

6.488E-11 

6.549E-11 

6.543E-11 

6.552E-11 

6.631E-11 

6.703E-11 

6.787E-11 

6.823E-11 

6.849E-11 

6.867E-11 

6.866E-11 

6.864E-11 

6.879E-11 

6.889E-11 

6.885E-11 

6.887E-11 

6.889E-11 

6.889E-11 

6.890E-11 

6.890E-11 

6 .889E-11 

6.890E-11 


-2.076E-14 
-2.086E-14 
-2.117E-14 
-2.128E-14 
-2 . 130E-14 
-2.149E-14 
-2.165E-14 
-2.181E-14 
-2.201E-14 
-2.220E-14 
-2.234E-14 
-2.242E-14 
-2.251E-14 
-2.275E-14 
-2.282E-14 
-2.309E-14 
-2.320E-14 
-2.343E-14 
-2.340E-14 
-2.342E-14 
-2.372E-14 
-2.400E-14 
-2.434E-14 
-2.449E-14 
-2.461E-14 
-2.469E-14 
-2.468E-14 
-2.467E-14 
-2.473E-14 
-2.478E-14 
-2.476E-14 
-2.477E-14 
-2.478E-14 
-2.478E-14 
-2.478E-14 
-2.478E-14 
-2.478E-14 
-2.478E-14 
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Appendix  B 
Table  B-II 


b  (E  (kev)) 
n  n 

Fit  Constants  for  Electron  Energy  Spectra 
from  Arbitrary  Neutron  Energies 

bo  bi 

b2 

bo 

6.093E-04 

1.076E-08 

-3.036E-13 

bl 

1.799E-05 

-2.364E-10 

5.214E-15 

b2 

-6.244E-08 

8.373E-13 

-1.903E-17 

b3 

6.847E-11 

-9.915E-16 

2.415E-20 

b4 

-2 .460E-14 

3.797E-19 

-9 .962E-24 

Appendix  C 


Electron  Angular  Distribution  Dependence 
on  Electron  and  Neutron  Energies 

To  show  the  dependence  of  the  electron  angular  distribution  on 
both  the  electron  and  neutron  energies,  this  appendix  presents  a 
series  of  plots.  The  data  for  each  plot  are  the  results  of  the 
monte-carlo  simulation  described  in  Chapter  III.  Each  plot  presents 
the  0-50  keV  electron  energy  group-0  and  the  500-550  keV  electron 
energy  group  group-0  for  twelve  neutron  energy  groups.  The  scales 
are  linear  on  both  axes  and  present  the  probability  of  observing  an 
electron  created  by  the  appropriate  neutron  of  a  specific  energy  as 
a  function  of  angle.  The  angle  is  defined  in  radians  from  0-tt  where 
a  0  angle  indicates  scatter  in  the  direction  of  the  neutron  velocity 
vector  and  a  ir  angle  indicates  backscatter  against  the  neutron 
velocity  vector. 

Isotropic  decay  in  the  center  of  mass  reference  frame  was 
assumed,  so  the  low  energy  neutrons  should,  and  do,  produce  such  an 
isotropic  distribution.  The  higher  energy  neutrons  should  skew  the 
faster  electrons  more  forward,  while  the  slower  electrons  remain 
nearly  isotropic. 

The  reason  for  this  co-dependence  is  that  the  resulting  electron 
angular  distribution  depends  on  the  velocity  addition  of  the  neutron 
and  the  electron  velocities.  Since  the  center  of  mass  energy 
distribution  remains  constant  and  the  center  of  mass  decays 
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isotropically,  the  lab  frame  will  observe  skewed  angular  distributions 
that  peak  in  the  direction  forward  to  the  perpendicular  of  the 
neutron  velocity.  The  energy  spectrum  in  the  lab  frame  will  also  be 
skewed  to  higher  energies  and  change  the  energy  distribution.  Since 
the  expected  number  of  electrons  of  a  specific  lab  frame  energy 
change,  so  will  the  overall  distribution  of  the  angles  those  elec¬ 
trons  will  be  seen.  Figures  C-l  through  C-12  show  the  numerical 
support  from  the  monte-carlo  simulation. 

Comparing  the  angular  distributions  of  all  electron  energy 
groups  caused  by  a  specific  neutron  energy  with  each  other,  it  is 
obvious  that  the  angular  distributions  are  different  except  for  those 
distributions  from  the  very  low  energy  neutrons.  Looking  at  the 
angular  distribution  of  a  specific  electron  energy  group  as  the 
neutron  energy  changes,  it  is  observed  that  this  distribution  also 
changes.  The  higher  energy  electrons  are  skewed  more  forward  with 
increasing  electron  energy,  while  the  lower  energy  electron  angular 
distributions  are  squashed,  but  remains  fairly  isotropic.  These 
observations  indicate  that  the  electron  angular  distribution  from 
neutron  decay  depends  on  both  the  electron  and  the  neutron  energy. 


PCE) 

0.000  0.012  0.024  0.036 


0.0000  0.78S4  1.5708  2.3562 

PHI  (RHDIRNS) 


3.1416 


Fig.  C-l  0-50  and  500-550  keV  Electron  Group  Angular 
Distributions  from  1.125E-03  keV  Neutron 


PHI  (RADIANS) 

Fig.  C-2  0-50  and  500-550  keV  Electron  Group  Angular 

Distributions  from  5.830E-01  keV  Neutron 


3 


0.0000  0.7854  1.5708  2.3562  3.1416 

PHI  (RADIANS) 

Fig.  C-3  0-50  and  500-550  keV  Electron  Group  Angular 

Distributions  from  2.478E+01  keV  Neutron 


0.0000  0.7854  1.5708  2.3562  3.1416 


PHI  (RADIANS) 

Fig.  C-4  0-50  and  500-550  keV  Electron' Group  Angular 

Distributions  from  5.503E+02  keV  Neutron 
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PHI  (RADIANS) 


Fig.  C-6  0-50  and  500-550  keV  Electron  Group  Angular 

Distributions  from  4.066E+03  keV  Neutron 
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Fig.  C-8  0-50  and  500-550  keV  Electron  Group  Angular 

Distributions  from  8.187E+Q3  keV  Neutron 


i 


P(E) 


3.1416 


0. 0000 


0.7854  1.5708  2.3562 

PHI  ( RADIANS) 


Fig.  C-9  0-50  and  500-550  Electron  Group  Angular 

Distributions  from  1.000E+Q4  keV  Neutron 


PHI  (RADIANS) 

Fig.  C-10  0-50  and  500-550  Electron  Group  Angular 

Distributions  from  1.282E+04  keV  Neutron 
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0.024  0.036 


3.1416 


0.0000  0.7854  1.5708  2.3562 

PHI  (RADIANS) 

Fig.  C-ll  0-50  and  500-550  keV  Electron  Group  Angular 
Distributions  from  1.492E+04  keV  Neutron 


0.0000  0.7854  1.5708  2.3562  3.1416 

PHI  (RADIANS) 

Fig.  C-12  0-50  and  500-550  keV  Electron  Group  Angular 

Distributions  from  1.964E+04  keV  Neutron 


Appendix  D 

Angular  Electron  Density  Along  a  Magnetic  Field  Line 


Figure  D-l  presents  the  view  of  an  arbitrary  neutron  decaying 
on  a  magnetic  field  line.  \J)  is  the  angle  of  incidence  of  the  neutron 
on  the  field  line.  0  is  the  angle  between  the  neutron  velocity 
vector  and  the  electron  velocity  vector.  And  a  is  the  angle  between 
the  electron  velocity  vector  and  the  magnetic  field  line.  The 
electron  distribution  as  a  function  of  0  is  known,  however,  the 
electron  distribution  as  a  function  of  ot  must  be  found  to  provide 
details  of  the  trapped  distribution  of  electrons. 


Fig.  D-l  Neutron  Decay  Geometry  Around  a  Magnetic  Field  Line 
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Defining  a  as  the  angle  from  the  direction  of  the  magnetic 
field  line,  and  $  as  the  angle  around  the  base  of  this  vector,  two 
spherical  coordinate  frames  are  defined.  To  calculate  the  distribu¬ 
tion  of  electrons  as  a  function  of  6,  a  coordinate  transformation 
is  required.  Figure  D-2  shows  the  relationship  between  the  spherical 
coordinate  systems.  The  equations  expressing  the  coordinates  in 
each  frame  are: 

x^  =  r  cosip  sin0  x^  =  p  cosg  sina 

=  r  simp  sin0  x^  ■  p  sing  sina 

x^  =  r  cos0  x^  =  p  cosa 


/ 

*, 

X, 


Fig.  D-2  Geometry  Axes 


The  transformation  between  these  coordinate  systems  is: 


X£  *  x^  cosip  +  x^  simp 
x3  *  -x£  simp  +  Xj  cos\p 
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These  relations  result  in  the  relation  that. 


cos0  -  -  sin$  sina  sinip  +  cosa  cosip  (D-l) 

which,  when  solved  for  6  equals, 

0  =  Arccos  (cosa  cosip  -  sin0  sina  sinip)  (D-2) 

For  a  given  a  and  ip,  the  values  of  0  can  be  found  by  iterating 
0  in  equation  D-2  from  0-2tt.  Summing  the  contributions  of  each 
electron  distribution  as  a  function  of  0,  each  increment  can  be  used 
to  calculate  the  electron  density  at  that  point.  The  result  is  the 
distribution  of  electrons  at  an  angle  a  to  the  magnetic  field  line. 

In  this  case,  where  N(0)  is  a  polynomial  distribution  of  0,  the 
number  density  in  N(0)  must  be  accounted  for  by  assuming  that  the 
electrons  are  isotropic  in  (p  and  by  dividing  N(0)  by  2nsin0.'  This 
calculation  determines  the  number  of  electrons  in  an  area  d0  by  dip. 
Since  no  simple  closed  form  solution  of  the  integrals  required  to 
solve  this  problem  exist,  the  solution  must  be  found  numerically  as 
described  in  the  previous  paragraph. 
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A  A 


Appendix  E 


User's  Manual  and  Subroutines  for 
Electron  Flux  Calculations  from  Prompt  Neutron  Spectra 

This  appendix  presents  the  subroutines  added  to  AFWL's  PROMPT 
radiation  code.  The  subroutine  DECAY  calculates  the  electron 
spectrum  and  flux  from  neutron  decay.  DECAY  uses  subroutine  EMAX  to 
calculate  the  maximum  electron  energy  from  a  specific  energy  neutron. 
Subroutine  RITDK  prints  the  electron  spectrum  and  flux  on  the  output 
file.  Subroutine  ANGLE  calculates  the  electron  angular  distribution 
as  a  function  of  electron  energy  and  angle  from  the  neutron  velocity. 

DECAY  requires  a  38  group  neutron  energy  spectrum,  array  SRSN, 
as  input.  It  defines  its  own  50  keV  electron  energy  group  spectrum 
in  array  EBINN  and  calculates  the  electron  energy  spectrum  from  the 
decayed  neutrons  in  array  SRSE.  The  variable  EFLUX  carries  the  value 
of  the  total  electron  flux.  Besides  the  neutron  source  spectrum,  no 
other  information  is  necessary.  A  listing  of  the  subroutines  is 
provided  on  the  following  pages. 

Subroutine  ANGLE  requires  the  same  38  group  neutron  energy 
spectrum  and  50  keV  electron  group  as  DECAY.  The  volume  of  data  is 
large  and  is  written  on  a  dummy  tape  called  TAPE9.  A  copy  of  this 
subroutine  is  provided  after  the  previous  mentioned  subroutines. 
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